UV-controlled physical and chemical structure 
of protoplanetary disks 
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Abstract We study details of the UV radiation trans- 
fer in a protoplanetary disk, paying attention to the 
influence of dust growth and sedimentation on the disk 
density and temperature. Also, we show how the dust 
evolution affects photoreaction rates of key molecules, 
like CN and CS. 
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1 Introduction 

Ultraviolet (UV) radiation is an important factor in 
the physical and chemical evolution of protoplanetary 
disks. It heats up a disk atmosphere, affecting both 
its structure and an emergent spectral energy distri- 
bution, and also controls the molecular content of the 
upper disk. A UV part of the spectrum comes primar- 
ily not from the star itself but from the inner region of 
an accretion flow (accretion shock, accretion column, 
etc.). As a result, the shape of the UV continuum and 
UV emission lines vary from a star to a star. Unlike 
the black body (BB) part of the stellar spectrum, its 
near UV part may significantly depend on details of a 



In this paper, we consider possible effects of a UV 
irradiation on the structure of a protoplanetary disk, 
in terms of density, temperature, and molecular com- 
position. This influence is related to disk heating by 
the central source and to the rates of photorcactions. 
As both factors depend on the overall distribution of 
the UV field in the disk, it is reasonable to expect that 
they are most influenced by the grain evolution (growth 
and sedimentation) which marks the very initial stage 
of planet formation. 

Radiative transfer (RT) modeling of a protoplane- 
tary disk is still a challenging problem (due to high op- 
tical depths and large variations of physical conditions), 
and special effort s are being made to d evelop appropri- 
ate methods fsee IPascucci et ah ( 2004 ) for benchmark- 
ing of RT codes). A large number of protoplanetary 
disk models has been created over the last two decades. 
Some of these models are focused on detailed (2D/3D) 
treatment of RT in dust continuum in order to describe 
the disk therm al structure a nd its spectral ene r gy dis - 



tribution (e.g. IWold (120031) . iNomura fc Millail (|2005l ) 



Dullemond fc Dominikl "( 20041 )). The goal of other mod 
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observations of T Tau stars (jAvresll201 



els is to study in detail the chemical and micro-physical 
struct ur e of protoplane t ary d is ks (e.g. pemenqv et al 
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(120071) . 

These simulations provide physical quantities (ioniza- 
tion degree, separate dust and gas temperatures) which 
control the dynamics of protoplanetary disks. These 
simulations are also extremely important to interpret 
existing and future observations of dust and molecular 
emission. 

Our goal is to develop a "balanced" model of a pro- 
toplanetary disk which would be relatively simple, but 
at the same time powerful enough to allow direct ob- 
servational verification. The model will be balanced in 
terms of complexity /reliability between RT and micro- 
physics treatment. With this model we plan to study 
the influence of different physical processes on the phys- 
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ical/cliemical structure and observational appearance 
of a protoplanetary disk. In the current version of 
the model we restrict ourselves to ID vertical RT but 
wit h scattering and deta iled frequency dependence (un- 
like IWoitke et all |2009l) where 2D RT is adopted but 
with only a few freque ncy channels). As shown in 
Dullemond et al.l (|2002l ). good frequency coverage is 
more important than angular ray coverage, that sup- 
ports our approximation. 



2 Disk Model 

To find mutually consistent density and temperature 
distributions in a disk one needs an iterative process. 
As a starting point, we use gas and dust distributions 
from IWiebe et al.l (|2008l ) . With these distributions be- 
ing set, we can solve a RT problem and find thermal 
structure of the disk, accounting for the external radia- 
tion (a central star with UV excess and diffuse galactic 
UV background). The disk temperature distribution is 
found by means of a two-stream model of the RT with 
detailed frequency grid (~ 200 frequencies covering a 
range from 100 nm to 1 mm). 

We assume that gas and dust are thermally coupled, 
while the temperature is determined by the dust ther- 
mal emission, absorption and isotropic scattering. The 
radiation field is separated into two streams, going from 
the surface of the disk to the midplane and in the op- 
posite direction. Such simplified consideration leads to 
the following RT equation in the z-direction: 
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^] is the mean intensity. 



where [ erg • cm ^ • ; 
X;y[cm^^] is the extinction coefficient, p is equal to 1/3 
or 1/4 for the cases of Eddington and Schwarzschild- 
Shuster approximations, respectively, and S^, is the 
source function: 
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(2) 
and 

are the absorption and scattering coefficients, 
Xv- The RT equation has to be solved si- 



Here B^(T{z)) is the Planck function, K^[cnr 
cr^ [ cm~ 

multaneously with the energy balance equation which is 
written under the assumption of radiative equilibrium: 
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The RT equation is solved using the analog of the 
Feautrier method with the triagonal (Thomas) algo- 
rithm for a hypermatrix. To carry out the iterations 



between mean intensity and the temperature T we 
employ the linear approximation of the Planck func- 
tion B^{T). The current approximation of T{z) allows 
to find the updated mean intensity distribution Ju{z). 
The solution of the RT equation in z-direction is re- 
peated for a number of radial distances i?, which gives 
the complete temperature distribution T{R,z). The 
dust opacities are calc ulated from the Mie theory (a s- 
trosilicate dust graines, Weingartner fc Drainj ( 2001 )). 

Once the thermal structure T{R, z) of the disk is 
known, we find the gas density structure p{z) at each 
R by integrating the equation of vertical hydrostatic 
equilibrium 
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where P is the gas pressure, fj, is the mean molecular 
weight (/^ = 2.3 for H2 and He mixture), nip is the pro- 
ton mass. Iterations are used to obtain self-consistent 
density and thermal distributions. In order to com- 
pletely solve the problem we need to specify a surfac e 
density S(i?) which was taken from lWiebe et al" ( 20081 ). 

The code has been extensively tested, first, for cases 
allowing analytical solutions. We considered a) an opti- 
cally thin media with the only heating source (radiation 
or accretion), b) a media with arbitrary optical depth 
irradiated by blackbody radiation field, and c) one- 
frequency case for a media with arbitrary optical depth, 
accretion heating and incident fiux. S econd, our model 



has be en compared with the model bv lDullemond et al 
( 2002r) . In their model, the direct stellar irradiation is 



considered to be the separate dust heating source, and 
the solution of the RT equation is being sought for the 
diffuse field. Our model produces colder upper layers, 
while midplane temperature is somewhat higher. We 
interpret the difference in the upper layer temperature 
for the two models as a possible consequence of under- 
estimated direct stellar irradiation in our model. At 
the same time, higher midplane temperatures in our 
model seem to be more realistic due to more accurate 
computation of dust radiative heating. 

We consider two dust models and four representa- 
tions for the spectrum of the central object. In model 
A5 dust is assumed to be well mixed with gas, with the 
mass ratio of 0.01. In model GS (growth and sedimen- 
tation) the dust distribution differs significantly from 
that of gas because of grain growth and sed imenta- 
tion iBirnstiel et al.ll201ol : IVasvunin et al]l201ll ). While 
the average grain size is 10~^ cm in model A5, in the 
midplane it increases because of dust evolution up to 
3 ■ 10^^ cm at the disk periphery (i? ~ 500 AU) and 
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10~^ cm in tlic inner disk region (we only consider lo- 
cations with R > 1 AU) . Large grains sediment to the 
disk midplane, which causes spatial variations in the 
dust-to-gas mass ratio. The ratio grows to few times 
lO^'^ in the midplane and decreases to 10~^ in the disk 
atmosphere. 

The stellar spectrum is assumed to be that of a black 
body with = 4000 K. We add a UV excess to this 
spec trum, describ ed by the scaled interstellar UV field 



Draine 



19781). a smoo thed BP Tau spectrum (JB; 



Kravtsova fc LamzinlboOS ). and a black body spectrum 
with temperature of 20000 K (J20). In all cases ex- 
cess UV continuum is added only for A < 4000 A. All 
added UV contributions are normalized to have the 
same mean intensity at A = 4000 A, as the star spec- 
trum. Also, we consider a case without a UV excess 
( J4) . Model designations are listed in Table [TJ 



3 UV irradiation and disk structure 

The influence of grain growth and settling on the disk 
structure has been studied in the literature a number 
of times 



(e.g. 



D'Alessio et al. 2006: Aikawa fc Nomura 



20061) . The general conclusion of these studies is that 
dust growth causes lower temperatures in the upper 
disk. The stellar radiation is mostly absorbed by small 
grains, and their number diminishes due to coagulation. 
On the other hand, dust sedimentation makes the up- 
per disk more transparent and hotter. The resultant 
temperature is thus defined by the net effect of dust 
evolution. 

Another factor which is important for disk heating 
is the UV continuum. In Figure [1] we show vertical 
temperature profiles at distances of 1 AU, 33 AU, and 
508 AU from the central source for models GS (left col- 
umn) and A5 (right column). For simplicity no itera- 
tions for the density structure has been performed for 
data on this figure. These profiles are quite different 
for the cases of pristine and evolved dust. In model 



Table 1 Model designations. 



Dust models 



GS Growth and sedimentation 
A5 Unevolving dust properties with 
an average size of 10 cm 
Incident UV spectra 
"J4 Black body spectrum with T, = 4000 K 
J20 Black body spectrum with = 20000 K 
JB Scaled BP Tau spcctum 
JD Scaled interstellar Draine UV field 



A5, where dust has the standard ISM properties, ver- 
tical temperature profile only weakly depends on the 
shape of the UV spectrum of the central source, while in 
model GS the dependence is more significant. If no UV 
continuum is taken into account, the disk with evolved 
dust (GS-J4) is indeed somewhat colder than the disk 
with pristine dust (A5-J4). However, if there is UV 
continuum, the disk with evolved dust is hotter for con- 
sidered cases. Temperature difference between various 
spectra is 120 K in model GS at 1 AU and less than 60 K 
in model A5. A similar trend is observed at other radii 
as well. Deeper in the disk, temperature is less sensitive 
to the dust model. Quite expcctedly, midplane temper- 
atures are somewhat higher in the more transparent 
disk with evolved dust. It is seen that the temperature 
range is bracketed by J4 and J20 spectra so further we 
consider only these cases. 

In Figure [2] we present results for the disk in hy- 
drostatic equilibrium, illuminated by the central object 
with spectra J4 and J20. In addition to models GS and 
A5, we also show density and temperature structure for 
the model with the same size distribution as in model 
A5 but with the upper size limit of 1 mm. It is seen that 
the upper limit of the grain size distribution does not 
change the density profile significantly, and correspond- 
ing curves in Figure [5] look very similar over the density 
drop of seven orders of magnitude. Model with evolved 
dust and without UV continuum (GS-J4) follows this 
trend as well. However, the curve for model GS-J20 
goes above all the other curves, indicating that the disk 
is puffed up in this model. It should be noted that 
curves for model GS end at lower heights than other 
curves. Dust density is very low at greater heights, 
and assumptions of our model break down. Obviously, 
a more detailed consideration is needed with separate 
treatment of dust and gas temperatures. 

Temperature profiles are more diverse. Models with 
an upper grain size limit of 1 mm are the coldest ones, 
in agreement with previous studies. Again, the im- 
portance of UV continuum is seen. For example, at 
100 AU it makes disk warmer by 5 K in 1 mm model, 
by 10 K in model A5, and by ~ 20K in model GS. 
This difference may not be strong enough to affect disk 
chemical composition, but it is of crucial importance 
for line transfer modeling. Note also that temperature 
gradient in model with evolved dust is much greater 
than in well-mixed models. This may cause differences 
not only in line intensities but also in the underlying 
disk molecular content. 



4 UV irradiation and disk chemistry 

UV irradiation of the disk coupled to the dust evolu- 
tion also should have a profound effect on its molecu- 
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Fig. 1 Vertical temperature distributions at 1.4 AU (top row), 33 AU (middle row), and 508 AU (bottom row) for evolved 
dust (left column) and pristine dust (right column). 



lar structure. To illustrate this, we use the disk mod- 
els described above to compute photodissociation rates 
for CS and CN molecules (Figure [3|), for which de- 
tailed frequency-dependen t cross-sections arc available 
( van Dishoeck et al hood) . Again, we only show results 
for J4 and J20 spectra. 

Obviously, photodissociation is much less effective 
in models without UV excess. Photodissociation rate 
of CS is lower by six orders of magnitude, while pho- 
todissociation rate of CN is lower by nine orders of 
magnitude. Also, these rates significantly depend on 
dust parameters, but this dependence is different for 
different molecules. For example, CN photodissocia- 
tion rate reaches value of, say, 10~^^ s~^ at z w 30 AU 
in model A5- J20 and at z « 20 AU in model GS-J20. 
For CS photodissociation rate the corresponding height 
range is almost twice as large. Such a different be- 
havior is related to wavelength dependence of reaction 



cross-sections. Molecules of CN are dissociated mostly 
by photons with A ^ 1000 A, while CS molecules can 
be destroyed by photons with large A ^ 1500 A. In 
model GS emission with longer wavelengths penetrates 
deeper into the disk due to dust settling, while emission 
with shorter wavelengths is effectively absorbed both in 
model A5 and in model GS. 

5 Conclusions 

In this paper we study the influence of UV continuum 
on the physical and chemical structure of protoplan- 
etary disks. Disk parameters (thermal structure and 
chemical photoreaction rates) may sensitively depend 
both on the mere presence of UV continuum, and on its 
exact shape. Thus, the UV radiation affects disk ob- 
servational appearance both in terms of continuum ob- 
servations and molecular line observations. Influence of 





Fig. 2 Vertical density and temperature distributions at 10 AU (left column) and 100 AU (right column) for various dust 
models and stellar spectra. 
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Fig. 3 Vertical profiles of CN and CS photodissociation rates at 100 AU for various dust models and stellar spectra. 



UV irradiation onto the disk structure gets stronger as 
dust particles grow bigger at the initial phase of planet 
formation. Thus, in order to interpret observations of 
a particular star+disk system its UV spectroscopy in 



the range of 1000 — 3000 A is greatly important, which 
makes future mission s lik e WSO-UV highly desirable 
(jShustov et al.ii200a . i201lh . 
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